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Abstract 



We have developed a method for calculating the two-point correlation function of non- 
linearly evolved mass and collapsed halos in the Press-Schechter formalism. The nonlin- 
ear gravitational interaction is treated as the sum of various individual spherical top-hat 
clustering. Because no collapsed halo of mass M can exist in initial regions (or top-hat 
spheres) of mass less than M, the bias that massive halos have stronger correlation than 
the background mass can be naturally introduced. 

We apply this method to derive constraints on popular dark-matter models from the 
spatial number density and the correlation function of C IV absorption systems in QSO 
spectra. Considering C IV systems should be hosted by collapsed halos, one can obtain 
an upper limit to the threshold mass of the collapsed halos by requiring their number 
density to be larger than that of observed C IV systems. On the other hand, in order to 
explain the observed clustering of C IV systems, a lower limit to the threshold mass will 
be set for the hosting halos. We found that the standard cold dark matter (SCDM) model 
and the low-density flat universe with a cosmological constant Aq (LCDM) are consistent 
with the abundance and clustering of C IV systems. However, the two cold-plus-hot 
dark matter models (CHDMs) with the cosmological parameters {Qc + ^b)/^h — 0.7/0.3 
and 0.8/0.2, respectively, have difficulty passing the two tests simultaneously. In these 
models, in order to have enough collapsed halos to host C IV systems, the threshold mass 
of the halos can not be greater than 10^^ Mq. But in order to agree with the two-point 
correlation function on the scales of Av ~ 300 — 1, 000 km/s, the threshold mass should 
be larger than 10^^ Mq. 

Subject headings: cosmology: theory — dark matter — quasars: absorption lines 
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1. Introduction 

It has been recognized that the abundance, i.e. the spatial number density, of moderate 
and high redshift objects can put promising constraints on models of structure formation 
in the universe. For instance, an N-body simulation has shown that the "standard" 
cold-plus-hot dark matter model (CHDM) of Qc/^h/^b — 0.6/0.3/0.1, here c, 6 and h 
denoting cold, baryonic and hot respectively, lacks the perturbations necessary to form 
clusters before z > 0.5, although it can produce a proper number of clusters at ^ = 
(Jing & Fang 1994). In contrast, the low density, flat cold-dark- matter model (LCDM) 
can produce enough clusters at both z > 0.5 and z — 0. Hence, the abundance of 2; > 0.5 
clusters is a useful quantity for discriminating between the CHDM model and the LCDM 
model. Another example is the damped Lya systems which are rich in QSO absorption 
spectra at 2; > 2. Both the Press- Schechter formahsm and N-body simulations indicated 
that the "standard" CHDM model has difficulty explaining these systems (Mo & Miralda- 
Escude 1994; Ma & Bcrtschinger 1994; Klypin et al. 1995). A simulation of the Lya forest 
has shown that the predicted number of forest lines is less than what is observed by at 
least one order of magnitude even if extreme parameters are used in this model (Bi, Ge 
& Fang 1995). 

Among the various samples of high-redshift objects, those selected from absorption 
systems in QSO spectra are relatively uniform and numerous and can be used to provide 
stringent constraints on different theoretical models. However, a common problem of 1-D 
samples is the lack of information on the size and geometry the objects. Calculating 3-D 
abundances from 1-D data requires assumptions about the objects' shape. Obviously, 
this leads to an uncertainty in the results. This difficulty is partly overcome in this study. 
We propose to test models not only by abundance, but also by clustering property of a 
high-redshift sample. Like abundance, clustering of collapsed halos depends also on their 
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mass and/or size, so adding a clustering test will reduce the unknown parameters. In 
some cases, conclusions could even be completely free from these parameters. 

The present analytical study is based on the Press- Schechter formalism (Press & 
Schechter 1974, hereafter PS) that has been found to be successful in describing the 
mean number density of collapsed virial halos with different mass threshold at high red- 
shifts. The spatial range in our study covers from redshift 2 to 4, and the correlation 
scales are from 300 to tens thousand km s~^, which are otherwise difficult to handle by 
current gas simulations. 

For our goal, we need a method to calculate the correlation function of collapsed halos 
in the PS formahsm. The existing PS theory does not yet tell us about the correlation 
function. In the first part of this paper, we will discuss how to extend the PS formalism 
to obtain the spatial correlation function of nonlinearly evolved mass and collapsed halos. 
A similar problem has been studied recently by Mo & White (1995). The idea is that 
the gravitational clustering from given initial density fluctuations can be approximately 
treated as many individual top-hat-evolved spherical regions. The mass correlation can 
then be deduced from the mass distribution within such regions. The bias of halo au- 
tocorrelations with respect to their masses is derived by assuming that no halo of mass 
M can exist in uncoUapsed regions of mass less than M. This approximation is found to 
be in good agreement with the linear approximation on large scales, and to be consistent 
with the empirical formalism on scales where the non-linear effects are significant (e.g. 
Hamilton et al. 1991). 

In the second part, we apply the developed method to C IV absorption systems in QSO 
spectra. We compare the abundance and clustering of collapsed halos in theoretical models 
with those given by real C IV observations. Among high-redshift absorption samples, only 
metal absorption systems and Lya forest lines can provide the statistics of both linear 
number densities and spatial correlation functions (Sargent, Boksenberg & Steidel 1988). 
Lya absorption lines are most likely from clouds with large size and low column density 
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which are neither viriahzed nor completely confined (Bechtold et al. 1994; Dinshaw et al. 
1994). Hence, we should not simply identify their hosts as collapsed halos. In contrast, 
metal systems are generally beheved to be from huge halos surrounding galaxies (Wolfe 
1993). Because stars were certainly formed in the systems, they must have undergone 
non-linear collapse. Therefore, metal absorption lines are probably the only available 1-D 
sample for our purpose. 

Observationally, metal absorption systems contain Mg II selected systems, C IV se- 
lected systems, and Lyman limit systems. However, the categories of absorbers identified 
by different systems are not orthogonal with each other (Wolfe 1993). The Lyman limit 
systems are almost indistinguishable from the Mg II systems, and most Mg II systems 
exhibit C IV, too. Most metal systems were obtained from ground observations. The C 
IV systems are usually found in the redshift interval 1.2 < z < 4.1, and the Mgll systems 
in 0.2 < ;s < 1.9. Because we are interested mainly in the physics at high redshift {z > 2), 
only C IV absorption systems will be investigated. 

It should be pointed out that, hke other studies of high-redshift objects based on the PS 
formalism, our goal is not to model the details of the C IV systems, instead it is to use the 
number of C IV systems as a lower limit to the number of corresponding collapsed halos. 
It is believed that the presence of C IV depends on the chemical abundance of heavy 
elements and on the ionization state of baryonic gas (Bergeron et al. 1994). Observations 
have shown an evolution in these chemical properties (Steidel 1990). It indicates that C 
IV systems should be located in areas in which chemical abundance evolution has already 
taken place. Therefore, C IV systems should be harbored in collapsed PS halos. Thus, 
a reasonable constraint on models of structure formation is that the predicted number of 
collapsed PS halos should be greater than that of observed C IV systems. 

In §2, we describe the PS method for calculating the spatial number density of collapsed 
halos, then we extend this idea to general uncoUapsed regions and derive an approximate 
expression for the correlation functions of mass and halos. In §3, the four SCDM, LCDM 
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and CHDM models of structure formation are tested based on their predictions of the 
abundance and clustering of collapsed halos. We discuss how observations can be com- 
pared with the predictions. Finally, §4 gives discussion and conclusions. 



2. Method 

2.1 Number density of collapsed halos 

The spatial number density of collapsed halos can be calculated from the standard 
PS theory. We define 5(x, z) to be the 3-D density fluctuation field of dark matter 
extrapolated to redshift z assuming linear evolution. A density field 5r(x), representing 
the smoothed fluctuation on scale R, can be derived from 5(x) by 



where the function W{R;xi — x) is the top-hat window for the comoving volume Vr — 
AttR^/3. In the = 1 Einstein-de Sitter universe, the variance of 5r evolves as cr^ oc 
(1 -I- ^)~^. The total mass within is M = VrPo, where po is the present cosmological 
density when the scaling factor of the universe is set to be unity at 2; = 0. 

For a Gaussian field, the fraction of the total mass po^x having fluctuations larger 
than a given 6c in an arbitrary spatial domain Ax is 



Therefore, if we take 5c — 1.686, the critical value for collapse in the top-hat evolution, 
Fr ■ poAyi should be identified as the sum of masses of all collapsed halos, each of which 
is greater than M — VrPq. The differential —^{FRpoAxjdM gives the total mass of 
collapsed halos in the range M to M -|- dM. Hence, if we define nc{M)dM to be the 
spatial number density of halos between M and M -\- dM, we have 




R 



(1) 




(2) 



d 



{FRpoAyi)dM = nc{M)dM ■ MAx, 



(3) 



dM 
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where we use the subscript c in ric to emphasize that it is for collapsed halos. 

Because the cloud-in-cloud problem has not been appropriately considered in Eq. (3), 
the original PS theory takes an ad hoc assumption that the above defined Uc should be 
multiphed by a factor of 2 (see discussions in Bond et al. 1991). The normalization 
/o°° 'nc{M)MdM = po can thus be fulfilled when (Ti^(i? — > 0) = oo. One has finally 

where erfc(x) is the complementary error function. The cumulative number density, 
iVc(M), of all halos greater than M should be 

POO 

N,{M) = / nc{M,)dMi. (5) 



IM 

The abundance of halos calculated from Eq. (5) has been verified by a number of N-body 
simulations (e.g. Efstathiou et al. 1988; Lacey & Cole 1994; Jing & Fang 1994). 

2.2 5c in the flat universe with a cosmological constant 

The critical value of collapse, Sc — 1.686, is derived in the Einstein-de Sitter model. 
Now, we calculate Sc in the flat A 7^ universe. In this case, the evolution of a spher- 
ical volume with mass M should be described by a Newtonian equation in the proper 
coordinate (Peebles 1984) : 
(Pr GM 1 

where A = (1 — f2)3i7Q. Eq.(6) can be integrated to give 

= 2GM{\ - + (1 - mlir' - ri), (7) 

where denotes the radius when the sphere reaches its maximum size. Prom Eq. (7), 
one has 

1 1 



t = r[2GM{ ) + {l-n)H',{r'-rl)]-'/'dr, 

Jo r Trr,. 



1 rrlvm 1 1 ._i/9 , 



(8) 



where cq = 2^7^/(1 - n)H^ri. 

If we define to be the cosmic time at which the sphere is at its maximum radius r^, 
the collapse time tc of this sphere is twice tm- Therefore, taking tm — tc/2 in Eq. (8), we 
find 



2 /"^ 1 1 
ic(co) = — / - y)-^co - - y)-^^^dy. 

Vl — "-no -^0 



(9) 



The redshift Zc that is a function of Cq can be found according to the z — t relationship in 
the fiat universe: 

2 



Hot 



-^sinh-h/i^(l + .)-3/^]. 



(10) 



3(1-0)1/2 Ly ^ 

At the time of recombination td, the radius of the sphere is much less than r^, so 
the integration in Eq. (8) can be approximated as 



td 



^|'ZA)3/2 _^ ^ + Co 2 ^Td ^5/2 



(11) 



a/I - ilifo 1-3 r^' ' 2co 5V^„ 
Again using the z — t relationship, one can write the density perturbation in the sphere 
at td, 5d = M/f r3(l + Zdfpo - 1, as 



(^d(co) = 



3co(l - n)H^r 



2^3 



l + Co3 



'i-n)co 



1/3 



1 + Zd 



(12) 



8G'7rr^(l + 2;d)3po cq 5 

Since depends only on cq, we can write the above cq as a function of Zd, i.e. cq = co{zd)- 
Using this initial density fiuctuation, one obtains the subsequent linear evolution at 
the epoch of collapse tc : 



5c{zc) = D{zc){l + Zd)5d = 



[Q(l + ..)^ + l-Q]V2^i±^ 

5 coizc) 



(13) 



'1 - Q)co{Zc 



nl/3 



(1 + z)dz 



[Q{1 + Zf + 1- 0]3/2 



where D[z) is the linear growing factor in the LCDM model. 
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Fig. 1 shows 6c{zc) vs. Zc for Q = 0.9, 0.7, 0.5, 0.3, 0.1. All of the Sc{zc) curves 
are approaching the traditional value 1.686; there is no noticeable difference among the 
thresholds at redshifts z > 1. This indicates that the trajectory of the top-hat collapse 
in the fiat A universe is very well described by simply the Einstein-de Sitter universe at 
high redshift, so we will take 5c = 1.686 in the following calculations. 

2.3 Number density of uncoUapsed regions 

As pointed out by Mo & White (1995), the nonlinear effect of gravitational clustering 
of dark matter can be treated as the sum of various individual top-hat spheres including 
both collapsed halos and uncoUapsed regions. An uncoUapsed region (sometimes called PS 
spherical region or uncoUapsed sphere; halo is usually a terminology for collapsed objects) 
corresponds to linear fluctuations not yet reaching the threshold 1.686 at the redshift 
considered. Therefore, wc can apply the same statistical method as that used by Press 
and Schechter for collapsed halos to uncoUapsed regions. Following §2.1, we first calculate 
the number density of uncoUapsed regions. 

The evolution of a spherical region of radius r = at recombination z = Zdis described 
by Eq. (6). In the Einstein-de Sitter case, the solution of Eq. (6) can be expressed 
analytically as 

for 5o > 0, and 

^ ^ _ cosh^), ^ = _^^i^(sinh^ _ ^)2/3 (15) 

10^0^ ' l + z 205o 

for do < 0, where ro = (1 -|- Zdjra and So = (1 -l- Zd)Sd are, respectively, the linear 

extrapolation of the radius and density contrast to 2; = 0. Note that the mass of this 

sphere is then given by M = ^(1 -|- Zd)^r^po(l + ^d)- 

Let's consider an arbitrary spherical volume of radius r at redshift z. Matter in this 

volume can come from various initial spheres of different ro and So- For a given r and 
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z, one can find the relationship between tq and 5q from Eqs. (14) and (15). This is 
^ — /(^-^7^)- Function f{x) is plotted in Fig. 2. The meaning of the figure is 
straightforward. When 5q — 0, one has r = + z), i.e. the evolution of this region 
is just comoving. An initial sphere with perturbation 5o > will evolve into the radius 
r < ro/(l + z) at rcdshift z. For 5q < 0, it will evolve into r > ro/(l + z). All initial 
spheres with 5q > 5c{l + z) result in r = 0, i.e. they are collapsed halos. 

One can consider the initial density field as a system consisting of many spheres, each 
of which has the same radius Tq but various density contrasts 5q. Similar to Eq. (2), the 
mass fraction of the spheres greater than 5q is given by 

/ ^ exp -— - (16) 

where (7^^ is the variance of density perturbation on scale tq. Because a spherical region 
of 6q and To will evolve into r at initial spherical regions oi 5 > 6q and Tq will evolve 
into radii less than r at z. Therefore, the comoving number density of the uncoUapsed PS 
spheres with mass in Mq Mq + dMo (or radii ro — > ro + tiro) at recombination, which 
are spheres with radius less than r at is given by 

5o 



n(Mo)dMo = -^l^erfc 



c^Mo, (17) 



.\/2a(ro) 

where we use the subscript in Mq to emphasis that it is for the uncollapsed regions, 
not only the collapsed halos discussed in §2.1. It should be pointed out that Eq. (16) 
contains spheres of both 5q < 5c{l + z) and 5q > 5c{l + z). Therefore, n(Mo)dMo, 
in fact, includes both uncollapsed regions and collapsed halos. Since initial spherical 
regions entirely cover the density field, n(Mo) should satisfy the normalization condition 
P — Io° '>T'{Mo)MQdMQ — po for all r and z. It is easy to verify this normalization from 
Eqs. (16) and (17). 

2.4 Mass-correlation function 

Since n{Mo)dMo gives the number density of uncollapsed spherical regions with mass 
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Mo — Mq + dMo and radii < r at 2;, the fraction of regions with radius r r + dr can be 
obtained by differentiating n{Mo) with respect to = Anr^/S. Defining m(Mo, V)dModV 
to be the number density of the regions with radius r^r + dr or V^V + dV, where 
dV — Airr^dr, we have 



m{Mo,V) = ^ri(Mo) 



Pol r 
-erfc 



So 



(18) 



Mo2dModV [V2a{ro)_ 

Therefore, the total number of such spherical regions in an arbitrary volume dVi should 

on average be given by m(Mo, V)dMQdVdVi. 

The mass correlation function ^(r) can be defined as the relative enhancement of mass 

density in the spherical shell r — > r + dr around dVi. Only the spheres with radius 

r — > r + can contribute to this enhancement. The mean enhancement of each Mq 

sphere is approximately described by its mass variance Mq(T^(Mo, 2;), where is the 

variance of the linear density contrast extrapolated to z. Therefore, the mass correlation 

function can be expressed as 

. . ^ /o°° m(Mo, V)dModVdVi ■ Mja^Mp, z) 
^^'''^^ PodV-podVi 

POO 

= / dMom{Mo,V)Vo^a^{Mo,z). (19) 
Jo 

The factors podV and podVi are the mean mass in the spherical shell r ^ r + dr and the 
volume dVi, respectively. Obviously, in deriving Eq. (19), we assumed that there is no 
correlation among the initial spheres, so the mass fluctuation in the shell is simply given 
by the sum of the individual components. 

We can check the approximation of Eq. (19) by comparing it with the linear approx- 
imation and other empirical non-linear formulae. First, because the mass correlation 
function is very well described by the linear approximation on large scales, Eq. (19) 
should be equal to the linear correlation function when r is large. The integrated mass 
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correlation function ^ that is defined by 



V Jo 



■V 



d_ 
dV 



(20) 



is shown in Fig. 3 for the SCDM, LCDM, CHDMl and CHDM2 models respectively 
dX z = 2.8. The solid lines in Fig. 3 represent our integrated correlation function from 
Eq.(19). The dotted lines in Fig. 3 are the linear approximation that matches Eq. (19) 
exactly on large scales. 

The dashed line at the upper left of Fig. 3, i.e. in the SCDM model, is the empirical 
correlation function of Hamilton et al. (1991) fitted to early N-body simulation results 
in the SCDM model. Our approximation has nonlinear clustering behavior similar to 
theirs except for higher correlations on small scales. Since there has been discussion as to 
whether early N-body simulations underestimated the correlation function on the smallest 
scales, a further judgment on these two approximations should be done by simulations 
with higher resolution. 

2.5 Correlation functions of collapsed halos 

We now consider the correlation function of collapsed halos. From Eq. (5), the number 
of collapsed halos in an arbitrary volume Vq is Nc{M)Vo. However, this is not the number 
within uncoUapsed spheres which we discussed in the last two sections, because it is 
impossible that a collapsed halo can form from uncoUapsed regions if the mass of the 
initial sphere is less than the mass of the collapsed halo. 

To avoid this difficulty, we require that inside an uncoUapsed sphere of mass Mq, 
the number density of collapsed halos is zero if their mass M is greater than Mq, or is 
proportional to Nc if their mass M is less than Mq. The number density of collapsed halos 
inside an PS sphere is thus given by 




for M < Mc 
for M > Mr 



0, 



(21) 
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where the constant A > 1 is introduced to maintain the normahzation condition 



Therefore, we have A = 2/erfc[5o(rc)/(7(rc)] and = (3M/47rpo)^/^- 

Obviously, Eq. (21) is a simphfied description of geometric bias (Kaiser 1984). Since 
the whole density field can be covered by spherical regions, every halo should be contained 
in a proper sphere. Massive collapsed halos can form only in massive uncoUapsed spheres. 
Therefore, the probability of finding a collapsed halo in a massive uncoUapsed sphere 
should be higher than in an average arbitrary volume by the factor A > 1. This effect 
will enhance the spatial correlation of massive halos. 

As in §2.4, we consider a typical spherical shell dV — Anr'^dr. The total number of 
uncoUapsed regions with radii r ^ r + dr and masses Mq Mq + dMo in the volume 
dVi is m(Mo, V)dModVdVi. The total number of collapsed halos greater than M in each 
Vq sphere is N'c{M)Vq, and the variance of the number is (A/'cH)^cr^(^o)- Therefore, the 
correlation function of collapsed halos with mass larger than M can be expressed as 



where we implicitly assume that in uncoUapsed regions, the collapsed halos have the same 
linear variance as the mass. 

3. Application to C IV systems 

We will concentrate on four models of structure formation: SCDM, LCDM, CHDMl 
and CHDM2. The parameters of these models are Usted in Table 1, including the Hubble 
constant h in the unit of 100 km s"^ Mpc~^, the cosmological density parameters Q, the 
cosmological constant Aq, and the quadrupole anisotropy Qrms of the cosmic microwave 
background radiation used for the normalization of the linear spectrum. The parameters 




(22) 



e(r;>M) 



dMomjMo, V)dVdV, ■ (A^Vp) V^(Mo, z) 
dVNc ■ dViNc 




(23) 
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are quite standard for these models. The hnear transfer functions of the SCDM and the 
LCDM models are taken from Bardeen et al. (1986) and that of the CHDMs from Klypin 
et al. (1993). For all four models, we have assumed the Harrison-Zel'dovich primordial 
power spectrum. It is worth pointing out that the Qrms normalization is compatible with 
the clustering of nearby galaxies in the LCDM, CHDMl and CHDM2 models, but is too 
high for the SCDM model. 

Part of the motivation for the two CHDM models comes from the possible non-zero 
rest mass of neutrinos. It has been claimed that the sum of the neutrino masses might 
be in the range 3.5 eV < Y^fT^m < 6 eV (e.g. Louis 1994). Therefore, the allowed range 
for the neutrino density would be 0.2 < Jl,^ < 0.3 when h — 0.5. The density parameters 
of the CHDMl and CHDM2 models are taken to be the maximum and minimum of this 
range. 

3.1 Abundance of possible hosts of C IV systems 

The redshift evolution of the number density, Nc{> M), of collapsed halos with mass 
greater than a given M are shown in Figs. 4a-d. The eight curves in each figure correspond 
to M = iQ^o+no.b ^j^j^ n = 0, 1...7 from top to bottom. 

As discussed in §1, QSO metal absorption systems should be hosted by collapsed halos, 
because only in such regions can the formation of stars and evolution of heavy elements 
be taking place. Therefore, the number of the observed C IV systems can be used as a 
lower limit to the number of collapsed holes. 

The number density of C IV systems is plotted as crosses in Fig. 4. The 3-D densities 
in the middle of the diagram were deduced from the observed line-of-sight number density 
(Steidel 1990) under the assumption that each absorber is spherical and has a radius of 
39 kpc (Crotts et al. 1994). To consider the effect of uncertainties in the radius, we 
also plot in Fig. 4 the 3-D number densities if the radius is 5 times less than, and 5 times 
greater than 39 kpc. These are the upper set of points (8.0 kpc) and the lower 
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set of points (195 kpc), respectively. The factor of 5 might be enough to account 
for uncertainties in the measurement of the radius. It would be safely to use the data 
corresponding to the 195 h~^kpc radius as the lower limit to the real number. 

Prom Fig. 4, one sees that the possible hosts of C IV systems in the SCDM and LCDM 
models should have a mass threshold of about 10^^'^ Mq. Say, that the number density 
of M > 10^^-^ Mq collapsed halos is great enough to fit with the observed C IV fines. 
Because the lower the mass threshold, the more the collapsed halos, any mass threshold 
less than 10^^-^ Mq is allowed by the observation. In the CHDMl and CHDM2 models, 
the number densities of collapsed halos are much less than those in the SCDM and LCDM 
models. In order to have enough collapsed halos hosting C IV systems, one has to choose 
much smaller mass thresholds. Figs. 4c and 4d showed that the assumed C IV halos 
should have a mass threshold as low as 10^^ Mq if the 195 kpc radius is used, or as 
low as 10^'^ Mq for the 39 kpc radius. This mass threshold is about one and a half 
orders of magnitude lower than that in the SCDM and LCDM models. We conclude that 
in the CHDM models, the halos that host C IV must have masses as, at least, low as 

10" Mq. 

3.2 Correlation functions of collapsed halos at high redshifts 

Using the approximation of Eq. (23), we present the correlation function ^{r, > M) for 
collapsed halos in the SCDM, LCDM, CHDMl and CHDM2 models in Fig. 5a-d. The 
mass threshold M is taken to be Mq. The eight curves in each model correspond 

to n = 0, ...7 from left to right, respectively. 

The histogram with error bars in Fig. 5 is the measured line-of-sight correlation func- 
tion of C IV absorption systems on scales from Af ~ 300 to 1,000 km s~^ (Sargent et al. 
1988). C IV pairs with velocity differences less than 300 km s~^ probably do not represent 
large-scale clustering, but the internal structures of the host galaxies. The peculiar ve- 
locities of the C IV absorbers may also significantly contaminate the spatial correlations 
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on scales less than the velocity dispersion. Therefore, the correlations below 300 km s~^ 
have not been shown in Fig. 5. Since the radii of the halos (39 kpc or even 195 
kpc), are much smaller than the scales being studied (> 300 km s~^), the 3-D correla- 
tion function is almost the same as that in 1-D. Therefore, one can directly compare the 
observed 1-D correlation length with theory. 

Fig. 5a shows that the correlation function of M > 10^^ M0 halos in the SCDM model 
provide a good fit to the observational data. The halos are also well within the upper 
hmit of the mass threshold, 10^^ Mq, derived in the last section. The same can be said 
for the LCDM model. In fact, SCDM and LCDM are indistinguishable in the top-hat 
evolution, because, as we argued before, a spherical mass has almost the same dynamic 
trajectory in all fiat universes with A < 0.8. Therefore, with the proper parameters, both 
SCDM and LCDM are consistent with the abundance and two-point correlation function 
of C IV systems. 

The two CHDMs are in trouble. According to the abundance fitting, the hosting 
halos of C IV systems should have masses as small as lO^^M©. However, the correlation 
function of such halos is much less than is observed. In order to match the observed 
correlation function, the halo masses must be at least lO^^M© (Figs. 5c and 5d). But the 
number of these halos is too few to account for the number of observed C IV systems. 
Therefore, there are no consistent parameters for CHDMl and CHDM2 under which both 
the abundance and correlation tests can passed. 

As pointed out by Heisler, Hogan & White (1989), the observed 1-D correlation function 
may not properly show the real-space correlation but may be amplified somehow by the 
peculiar velocity. According to their estimation based on a simple model, the amplitude 
of the correlation function on scales less than the peculiar velocity could be amplified by 
the factor of {ro/rd)^, where ro is the correlation length, 7 the correlation index and rd 
the cloud radius. Namely, the correlation length is amplified by the factor {ro/rd). In the 
CHDM models, the correlation length of IO^^Mq halos is less than 0.1 Mpc at redshift 
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2.8, so the amplification factor does not exceed 3. Tfiereforc, this mechanism is far from 
resolving the discrepancy which are found on scales of 3 - 5 h~^ Mpc in the CHDMs. 

4. Discussion and conclusions 

Using the approximate expressions of correlation function and halo's number density 
in the PS formalism, we have derived constraints on the masses of collapsed halos in 
four models of structure formation. Generally speaking, because the number density of 
collapsed halos is inversely dependent on the threshold mass above which the halos are 
selected, the abundance of C IV systems sets up an upper limit to the threshold mass 
for their hosts. On the other hand, the 2-point correlation function of collapsed halos is 
positively dependent on the threshold mass, so the C IV clustering observation sets up a 
lower limit to the threshold mass. Therefore, a model should be considered inconsistent 
if the lower limit is found to be larger than the upper limit. The two CHDM models are 
just such examples. Even when the size of the C IV absorbers is taken to be 5 times 
larger than that observed, the upper hmit (IO^^Mq) provided by the abundance is still 
inconsistent with the lower limit (lO^^M©) inferred from the correlation. 

We note from Fig. 3 that the mass correlation functions of CHDMs are much less 
than those of SCDM and LCDM. Yet, their corresponding halo correlations are greatly 
enhanced if we look at the same mass thresholds in Fig. 5. This is due to the large 
gravitational bias of massive halos in the CHDM models. Could it be possible that other 
biasing effects can further enhance the correlations, so the models can finally pass the 
correlation test? The answer probably is negative, because, expect gravitation, there 
are almost no available bias mechanisms for generating a higher correlation for less mass 
halos. Hydrodynamical processes generally are ineffective to produce any inhomogeneity 
on scales equal to or larger than 5 h~^Mpc, as the streaming velocity of gaseous component 
is only about a thousand km s~^ on average. 

Changing the density parameter of the hot dark matter from Qh — 0.3 (CHDMl) to 
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0.2 (CHDM2) does not relieve this difficulty much either. Though CHDM2 has a variance 
1.2 times that of CHDMl in the Unear evolution, it still cannot produce enough halos 
hosting C IV systems. 

Finally, as pointed out in §2.4, our approximation of the correlation function may give 
higher power on small scales than some N-body simulations; however, this does not affect 
the conclusion about the CHDMs models. If Eq. (23) really overestimates the correlation 
function on small scales, the difference between the models and the observation would be 
even more significant than is shown in Figs. 5c and 5d. In that case, the difficulty for the 
models would become more severe. 
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Table 1. Parameters of the cosmological models 





SCDM 


LCDM 


CHDMl 


CHDM2 
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0.5 


0.75 


0.5 


0.5 


^cb 


1 


0.3 


0.7 


0.8 










0.3 


0.2 


Ao 





0.7 










18.0 


18.0 


18.0 


18.0 
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Figure Captions 



Figure 1 Threshold 6c as a function of the collapsing redshift Zc for Q, = 0.9, 0.7, 0.5, 
0.3, 0.1, plotted from top to bottom. 

Figure 2 The function f{x) defined by 6o = Sc{l + z)f{r{l + z)/ro) according to Eqs. 
(14) and (15). 

Figure 3 The integrated mass correlation functions in the SCDM, LCDM, CHDMl and 
CHDM2 models. The solid hnes are calculated from Eqs.(19) and (20). The dotted 
hues are the hnear correlation functions, and the dashed hne in the SCDM model 
is from the empirical formula of Hamilton et al. (1991). 

Figure 4a The comoving abundance of collapsed halos with masses greater than M in 
the SCDM model. The eight curves correspond to M = lO^^+^'^^M©, n = 0, 1, .... 7, 
from top to bottom. The bold error bars in the middle are the number densities 
of QSO C IV absorption systems deduced from 1-D observational data and the 
measured radius 39 kpc. The data sets above and below the middle ones are 
the number densities when the radius is (39/5) kpc and (5 x 39) kpc, 
respectively. 

Figure 4b The same as Fig. 4a for the LCDM model. 
Figure 4c The same as Fig. 4a for the CHDMl model. 
Figure 4d The same as Fig. 4a for the CHDM2 model. 

Figure 5a The correlation function of collapsed halos with masses greater than M in 
the SCDM model. The eight curves correspond to M = lO^o+o-^" M©, n = 0, 1, 7, 
from left to right. The histogram with error bars is the measured correlation function 
of C IV systems. 
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Figure 5b The same as Fig. 
Figure 5c The same as Fig. 
Figure 5d The same as Fig. 



5a for the LCDM model. 
5a for the CHDMl model. 
5a for the CHDM2 model. 
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